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Abstract 

The sensitivity of the present CNGS beam to the sub-dominant <-> v e oscillations in the 
region indicated by the atmospheric neutrino experiments is investigated. In particular, we 
present a revised analysis of the OPERA detector and discuss the sensitivity to #13 of ICARUS 
and OPERA combined. We show that the CNGS beam optimised for v r appearance, will 
improve significantly (about a factor 5) the current limit of CHOOZ and explore most of the 
region sin 2 20i 3 ~ 0(1(U 2 ). 



1 Introduction 




There are currently intense efforts to confirm and extend the evidence for neutrino oscillations in the at- 
mospheric and solar sectors. In particular the disappearance of muon neutrinos originating from cosmic 
ray interactions provides a strong hint in favour of <-> v T oscillations. This indication will be 

tested in a conclusive manner at the next generation of long baseline experiments. In a simple two-flavour 
mixing scheme the atmospheric neutrino experiments constrain the mixing parameter in the region around 
(Am 2 , sin 2 20)=(2.5 x 1CT 3 , 1). Recently, the SNO experiment gave an unambiguous signature of the 
solar neutrino oscillations, excluding the sterile neutrino hypothesis 

In addition to the dominant <— ► v r component, sub-dominant oscillations involving also v e could 
be observed. In this three-flavour scenario the gauge eigenstates v a (a = e, /i, r) are related to the mass 
eigenstates v[ (i = 1, 2, 3) through the leptonic mixing matrix U a i- This matrix can be parametrized as: 

C12C13 S12C13 

U"(012,013,023,<5) = I -S12C23 - C12S13S23C" 5 Ci 2 C 23 ~ Sl2Sl3S23e i<5 Ci 3 S 2 3 1 (1) 

S12S23 - Ci 2 S 13 C 2 3e l5 -C12S23 - S 12 Si 3 C 23 e lS 

with S{j = sin Oij and Cij = cos Oij . 

Here the leading «-> v T transitions are driven by the angle sin 2 2023 an d the squared mass difference 
Am 2 3 . Sub-dominant <— > v e oscillations at the atmospheric scale (Am 2 3 ~ 2.5 X 10~ 3 ) are driven 
by the mixing angle $13. This angle is constrained by the CHOOZ experiment to be small (sin 2 2^13 < 

The next generation of long baseline experiments — MINOS f§, ICARUS § and OPERA @ — will 
search for a conclusive evidence of the oscillation mechanism at the atmospheric scale. In particular, the 
CNGS beam J^,^] will directly confirm the indication for neutrino oscillations through tau appearance. 

Although the precision that will be achieved on the parameters driving the dominant <— > v T after 
this first generation of long baseline experiments is well established Jll| , the sensitivity to the sub-dominant 

^ v e channel deserves further investigation, in particular for what concerns CNGS. This is due to the 
fact that only recently the baseline configuration of the CNGS detectors (ICARUS and OPERA Q) has 
been established. Moreover, a full <-> v e analysis with the OPERA experiment has not been carried out 
yet. 

In this paper we discuss the sensitivity of the ICARUS and OPERA experiments in searching for #13 
by using the present CNGS neutrino beam. The sensitivity of ICARUS to #13 has been evaluated in Q]. 
For OPERA, whose primary aim is the direct detection of v r appearance, a preliminary estimate of the 
sensitivity to y« *-* v e oscillations has been given in the 1999 Letter of Intent |ll| and is evaluated here in 
more details. The paper is organised as follows: in Section ^ we briefly describe the present CNGS neutrino 
beam, the ICARUS and OPERA experimental setups and their features; in Section || we discuss the ICARUS 
and OPERA sensitivity to #13 ; finally in Section || we draw our conclusions. 
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2 The present CNGS experimental program 



2.1 The CNGS neutrino beam 

In November 2000 a new version of the CNGS beam design was released @. In Tables g and [| the perfor- 
mance of the beam are given. 

Table 1: Nominal performance of the CNGS reference beam fl(i|| . 



Up (m 2 /pot) 


7.78 x 10~ 9 


CC events/pot /kton 


5.85 x 10~ 17 




0.8% 




2.1% 




0.07% 



Table 2: Number of v r charged-current interactions at Gran Sasso per kton and per year (shared mode) for 
different values of Am 2 and sin 2 (28) = 1 Q. 



Am 2 


v r CC interactions/kton/year 


1 x 10" 3 eV 2 


2.53 


3. x lO- 3 eV 2 


22.8 


5 x 10" 3 eV 2 


63.3 



The possibility of an increase of the neutrino beam intensity at moderate cost is being studied by the 
CNGS design group |Q and other groups working on the PS and SPS machines [ pS] . The SPS and its 
injectors have a crucial role for the LHC accelerator and considerable efforts and investment are going into 
the upgrade of these accelerators to meet the requirements of LHC. Most of the efforts are also directly 
beneficial to CNGS: for example, it is expected that the maximum intensity accelerated in the SPS will be 
7 X 10 13 rather than 4.5 X 10 13 pot/cycle Moreover, it has been recently pointed out that, taking 
into account that all the accelerators (LINAC, PSB, PS and SPS) are close to their limits, an increase in 
the proton beam intensity can be reached - pending further accelerator machine development studies - by 
reducing the beam losses between two sequential accelerators. These improvements could, at a moderate 
cost, increase the proton flux by as much as a factor of 1.5. 

Other scenarios for a further increase of the SPS extracted beam intensity has been outlined in pa], but 
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those are much more expensive and might turn out to be technically unrealistic. It is therefore prudent to 
state that an intensity increase for CNGS of a factor 1.5 is within reach. 

The CNGS energy spectrum has been optimized for v T appearance and the prompt v T contamination 
(mainly from D s decays) is negligible. On the other hand, in the CNGS beam, the expected v e contamination 
is relatively small compared to the dominant component, see Table |l| and allow to search for the sub- 
dominant oscillation *-* v e seeking after an excess of v e charged-current (CC) events. 

The systematic error associated with the v e contamination plays an important role for the Un <— > v e 
oscillation search, the statistical fluctuation of this v e component being the irreducible limiting factor. In 
turns, this uncertainty depends on the knowledge of the K yield. For sake of consistency with Refs. IO,E9[, we 



assume for this study a 5% systematic error on the overall u e flux. However, as we will show in Section 3.1 



due to the small number of selected events, especially in OPERA, a more conservative estimate of the 
systematics (up to 10%) would not change appreciably the experimental sensitivity. 



2.2 The ICARUS experiment 

The aim of the ICARUS Collaboration jl^] is the construction of a liquid argon time projection chamber 
(TPC) detector with a space resolution comparable to that of bubble chambers. 

In an ICARUS module a large liquid argon volume is contained in a cryostat equipped with wire readout 
planes and a high voltage plane at the center giving a long drift length. Ionisation electrons produced by 
charged particles crossing the liquid argon are collected at the readout planes. The arrival time of the 
different electrons allows the reconstruction of the particle track trajectory by the time projection technique. 
Each readout plane has three coordinates at 60° from each other and wire pitch of 3 mm. This corresponds 
to a bubble diameter of about 3 mm. 

The detector is a homogeneous tracking device capable of dE/dx measurement. Its high dE/dx 
resolution allows good momentum measurement and particle identification of soft particles. Electromag- 
netic and hadronic showers are, in fact, identified and fully sampled. This gives the energy resolution 
for electromagnetic showers of a(E)/E = 2>%/^jE{GeV) © 1% and for hadronic contained showers of 
a{E)/E = 16%/y 7 E{GeV) © 1%. ICARUS features excellent electron identification and e/7T° separation. 

The physics programme of the ICARUS experiment as described in the Proposal |l^] has been derived 
assuming a "baseline configuration" consisting of one T600 plus two T1200 modules for a total fiducial mass 
of 2.35 kton §. 



2.3 The OPERA experiment 

The OPERA experiment J§|;[lJ] is meant for a long baseline search for <-> v T oscillations in the CNGS 
beam fl4|| . The experiment uses nuclear emulsions as high resolution tracking devices for the direct detection 
of the r produced in the v T CC interactions with the target. Due to its excellent electron identification 
capability exploited for the reconstruction of the r — > e decay, OPERA is able to perform also a *-* v e 
oscillation search. 

OPERA is designed starting from the Emulsion Cloud Chamber (ECC) concept (see references quoted 
in H) which combines the high precision tracking capabilities of nuclear emulsions and the large mass 
achievable by employing metal plates as a target. The basic element of the ECC is the cell which is made 
of a 1 mm thick lead plate followed by a thin emulsion film. The film is made up of a pair of emulsion 
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layers 50 /im thick on either side of a 200 /im plastic base. Charged particles give two track segments in 
each emulsion film. The number of grain hits in 50 /im (15-20) ensures redundancy in the measurement of 
particle trajectories. 

For details on the event reconstruction both with the nuclear emulsions and the electronic detector, and 
the OPERA sensitivity to v T we refer to H|l|. The mass of the OPERA "baseline configuration" is 

1.77 kton fil| . In the following an effective mass of 1.65 kton has been assumed to account for the brick 
removal during data taking. 

Thanks to the dense ECC structure and the high granularity provided by the nuclear emulsions, the 
OPERA detector is also suited for electron and 7 detection ||. The resolution in measuring the energy of 
an electromagnetic shower in the energy range relevant to CNGS is approximately constant and is about 
20%. Furthermore, the nuclear emulsions are able to measure the number of grains associated to each track. 
This allows an excellent two tracks separation (better than 1 /im). Therefore, it is possible to disentangle 
single-electron tracks from tracks produced by electron pairs coming from 7 conversion in the lead. These 
features are particularly important for the Un <— ► v e analysis. 

The outstanding position resolution of nuclear emulsions can also be used to measure the angle of 
each charged track with an accuracy of about 1 mrad. This allows momentum measurement by using the 
Multiple Coulomb Scattering with a resolution of about 20% and the reconstruction of kinematical variables 
characterising the event (i.e. the missing transverse momentum at the interaction vertex p™ %ss and the 
transverse momentum of a track with respect to hadronic shower direction Qt)- 

3 Sensitivity to #13 
3.1 Event selection 

In this Section we discuss the event selection in OPERA and the expected number of signal and background 
events in both ICARUS and OPERA. 

The overall rate is defined as 

R(Amj 3 , sin 2 29 23 , sin 2 20 13 ) = S + + 5("mCC-u/ m NC) + ^(beam) + ^(NC) ; ^ 

where the signal, S, is the convolution of the flux [<j) v [E)) with the v e charged-current cross section 
(a^(E)), the U^ «-> v e oscillation probability {P u ^ Ue (E)) and the the signal efficiency (e s ^S na ^) 

S = Aj K (E) P u ^ Ue (E)ag c (E) e si § nal dE . (3) 

The normalisation A takes into account the effective target mass. The background coming from the 
dominant Un — ► v T oscillation channel where the produced r decays in one electron is 

B^ e) =Aj ^(E) P^ r (E)ag c (E) e^dE , (4) 

where e^ T ^ e ^ is the efficiency and is the v T charged-current cross-section. This background includes 
mainly taus decayed in the same lead plate as the primary vertex ("short decays") plus a small (4%) 
contribution from "long decays" where the kink has not been detected. 
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The background from events with the primary muon not identified is given by 

b(^CC-^NC) = A J ^{E)P u ^ i (E)a^{E)e^ CC ^ NC) dE , (5) 

where gC^CC— ►iyNC) j g ^ e probability for a z^CC to be identified as a f^NC and with another track 
mimicking an electron. 

The background coming from the v e beam contamination is 

5 (beam) = aJ <j) Ue (E) P Ve ^ Ue (E)vg C (E) e^^dE , (6) 
where g( beam ) is the efficiency for CC events from the v e beam contamination. 

The last background included, i?( NC ), from the decay of neutral pions created by neutral-current (NC) 
interactions, is 

=Aj ^(E) a* c (E) e^dE , (7) 

where e^ NC -* is the reduction efficiency. Further background from charge exchange process is negligible 

The OPERA analysis for <-> u e oscillation quoted in the 1999 Letter of Intent has been revised. 
The major improvements are the following: 

• the muon identification is improved from 80% to 95% [g. This means that the contribution to the 
background from f^CC events can drastically be reduced; 

• 7's produced at the primary vertex (7T° decays) and converting into the lead plate where the interaction 
occurred are identified if at least one of the below conditions is fulfilled 

1. the e^e - pair has an opening angle larger than 3 mrad. This cut is fixed by the angular resolution 
achievable by the emulsion films; 

2. if the e + e~ pair has an opening angle smaller than 3 mrad, we exploit the nuclear emulsion 
capability in measuring the number of grains associated with a reconstructed track. The rejection 
power of this cut has been computed by assuming for a minimum ionising particle 30 grains per 
emulsion sheet, while 60 grains are assumed for the e + e _ pair with opening angle smaller than 
3 mrad. A reconstructed track is classified as an e + e~ pair if the associated number of grains is 
larger than (30 + 3 X \/30). Conservatively we consider only the emulsion film downstream from 
the vertex lead plate. Note, however, that a significant improvement could be obtained including 
the grains in the subsequent emulsion sheets, crossed by the electrons before showering; 

• since neutrino oscillations create an excess of u e CC induced events at low neutrino energies, a cut on 
the visible energy is applied. 

Hence, the OPERA <-> v e dedicated analysis seeks for neutrino interactions with a candidate electron 
from the primary vertex with an energy greater than 1 GeV (this cut reduces the soft 7 component) and 
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a visible energy of the event smaller than 20 GeV (this cut reduces the prompt v e component of the back- 
ground). Moreover, we also apply a cut on the number of grains associated with the track of the candidate 
electron. The latter has a strong impact on the reduction of the background from f^CC and z^NC events 
and allowed for a softer cut on the electron energy. Finally, a cut on missing transverse momentum at the 
primary vertex pjp tss of the event is applied (p™ lss < 1.5 GeV) to further reduce NC contaminations and 
suppress r — > e background. 

For the brick finding, the vertex finding, the trigger efficiencies and the fiducial volume cut, we used 
the numbers quoted in the Status Report |Q . The product of these efficiencies is shown in Table || and is 
indicated with £. 

The final efficiencies (e) for both the signal and the background channels are shown in the second raw 
of Table || and the expected number of events is shown in Table ||. 

The search for z/„ <-> u e oscillations does not increase the emulsion scanning load by a large amount. 
Indeed, once the events have been located, the time needed to search for interesting topologies around the 
reconstructed vertex is small compared with the time needed to locate the events. On the other hand, for 
events identified as NC produced in the downstream half of the brick one should always remove the following 
brick to perform the electron identification. It implies that 13% more bricks have to be removed and scanned. 
Notice that in the preliminary estimate quoted in the 1999 Letter of Intent only half of the ECC mass has 
been used. 

Table 3: Signal and background efficiencies for the <-* v e oscillation search with OPERA. Notice that 
the t — > e efficiency includes the r — > evv branching ratio. 





^ e CC signal 


t — ► e 


^CC -> t^NC 


^NC 


f e CC beam 




0.53 


0.53 


0.52 


0.48 


0.53 


e 


0.31 


0.032 


0.34 x 10~ 4 


7.0 x 10~ 4 


0.082 



If we assume that both i>n <— ► v T an d fn <~ * v e oscillations occur simultaneously with oscillation pa- 
rameters A77I23 = 2-5 x 1CU 3 eV 2 , 023 = 45° and #13 G [3° 9°], the expected number of signal and 
background events both in ICARUS and OPERA are given in Table |]. The rates are normalised to 5 years 
data taking and assuming the nominal intensity of the CNGS beam. 

The ICARUS rates at CNGS from <-> v e oscillations and from background sources have been taken 
from Table 5.9 of Ref. [7J and scaled at the different oscillation parameters (Am| 2 , sin 2 2#23) sin 2 2^13) . 
Here a cut on the electron energy (E e > 1 GeV) and the visible energy (E v i s < 20 GeV) have been 
applied. In Table |4| the background from 7T° production in neutral-current events is negligible for ICARUS 
as assumed in Ref. M. This is not true for OPERA being the lead radiation length much shorter than 
the one of liquid argon. Nevertheless, this background can be further suppressed by a kincmatical analysis, 
see Section 3.2. On the other hand, the r — > e background in OPERA is suppressed by the capability of 
this detector to reconstruct exclusively the r decay topology. Finally, notice that in both experiments the 
dominant background remains the v e contamination of the beam. 
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Table 4: Expected number of signal and background events for both ICARUS and OPERA assuming 5 years 
data taking with the nominal CNGS beam and oscillation parameters Am^ — 2.5 x 10~ 3 eV 2 , 6*23 = 45° 
and i3 € [3°-=-9°]. 



013 


sin 2 20i3 


u e CC signal 


r — ► e 


z^CC -> i/„NC 


z^NC 


I'eCC beam 


ICARUS 














9° 


0.095 


27 


24 






50 


8° 


0.076 


21 


24 






50 


7° 


0.058 


16 


24 






50 


5° 


0.030 


8.4 


25 






50 


3° 


0.011 


3.1 


25 






50 


OPERA 














9° 


0.095 


9.3 


4.5 


1.0 


5.2 


18 


8° 


0.076 


7.4 


4.5 


1.0 


5.2 


18 


7° 


0.058 


5.8 


4.6 


1.0 


5.2 


18 


5° 


0.030 


3.0 


4.6 


1.0 


5.2 


18 


3° 


0.011 


1.2 


4.7 


1.0 


5.2 


18 



3.2 Kinematical analysis 

An increase of sensitivity to <— > v e oscillations can be obtained by fitting the kinematical distributions 
of the selected events. As an example the OPERA p™ tss distribution for signal and background channels is 
given in Fig. |l|. 

In order to disentangle <-> v e and *-* v T oscillations, the ICARUS Collaboration considered 
the transverse missing momentum. The Un <— > v e events are balanced, while <— > v T show a sizable 
transverse momentum due to the presence of two neutrinos in the final state. The visible energy and the 
missing transverse momentum are combined into a binned x^-?it, in which the three oscillation parameters 
are allowed to vary. 

By using the information contained in the Figs. 5.9 and 5.10 of Ref. JtJ, we have performed a binned x 2 ~ 
fit of these distributions and reproduced the results reported by the ICARUS Collaboration in Ref. j?J . The 
sensitivity at 90% C.L. to #13, corresponding to x 2 = Xmin + 4-6, has been evaluated under the assumption 
$23 = 45°. The exclusion plot is reported in Fig. |[ 

In the case of OPERA, the kinematical variables have been obtained starting from true quantities 
and applying a smearing according to the OPERA experimental resolutions, see Section |2.3| . By fitting 
simultaneously the E v i s , E e , p™ lss distributions , we obtained the exclusion plot at 90% C.L. shown in 
Fig. H under the assumption 023 = 45° . 

The use of differential distributions improves the sensitivity to sin 2 20 i3 of both ICARUS and OPERA 
by about 20%. The combined sensitivity at 90% C.L. of the ICARUS and OPERA experiments is obtained 
by minimising the x 2 defined as x 2 = xl( icarus ) + x\ opera) and is also shown in Fig. |. 
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Figure 1: p™ lss distribution for signal and background channels assuming #13 = 8°, Am| 3 = 
2.5 x 10- 3 eV 2 and 6*23 = 45°. 
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Figure 3: OPERA sensitivity to the parameter #13 at 90% C.L. in a three family mixing scenario, 
in presence of <-> v T with #23 = 45°. The sensitivity with the higher intensity CNGS beam is 
also given. 
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3.3 Results 



In the previous Sections we discussed the sensitivity achievable in searching for #13 with the present CNGS 
beam by combining the ICARUS and OPERA experiments. The improvements on sin 2 2#i3 over the CHOOZ 
limit in the region indicated by atmospheric neutrinos is about a factor 5 assuming 5 years data taking with 
the nominal CNGS beam. The limits at 90% C.L. on sin 2 29is and #13 for various experiments are given in 
Table |5| The limits achievable with an improved CNGS beam (1.5 times more proton) are also shown. 

It is worth noticing that the sensitivity of the combined experiments running with the present high 
energy CNGS beam is comparable to the one achievable with a dedicated low energy configuration running 
below the kinematical threshold for r production jH) and is obtained "for free" , as a by-product of the 
search for v T appearance. 

Clearly, we expect JHF phase 1 to improve significantly the sensitivity to sin 2 2#i3 with respect to the 
CNGS programme in case of a negative result. 

Table 5: Limits at 90% C.L. on sin 2 2#i3 and #13 for various experiments. 



Experiment 


sin^ 26»i3 


6*13 


CHOOZ 


< 0.14 


C 11° 


MINOS 


< 0.06 


< 7.1° 


ICARUS 


< 0.04 


< 5.8° 


OPERA 


< 0.06 


< 7.1° 


ICARUS and OPERA 






combined (CNGS) 


< 0.03 


< 5.0° 


ICARUS and OPERA 






combined (CNGSxl.5) 


< 0.025 


< 4.5° 


JHF [2(J 


< 0.006 


< 2.5° 



4 Conclusion 

The OPERA experiment was designed and optimised for the direct observation of the r decay topology 
and to search for v T appearance at the Am 2 atmospheric scale. Nevertheless, given the good electron 
identification, it turned out to be sensitive to the sub-dominant <— > v e oscillations in the region indicated 
by atmospheric results. Exploring the region of sin 2 2#i3 ~ 0(1O -2 ) is particularly interesting in connection 
with the discovery of CP violation in the leptonic sector pCy . Hence, an observation of sin 2 2#i3 in the region 
of 0(1CU 2 ) would place future CP violation searches like JHF phase 2 and the Neutrino Factories on firmer 
ground. We showed that by using the present CNGS beam and combining the ICARUS and OPERA 
experiments a sensitivity in sin 2 2^13 down to 0.03 (Am 2 3 = 2.5 X 1CU 3 eV 2 and 623 = 45°) can be 
achieved. It can be further improved if a higher intensity (6.8 X 10 19 pot/year) will be reached. 
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